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INTRODUCTION 


This  report  documents  the  analytical  study  of  the  temperature  changes  to  the  propellant  of 
a  tank  round  after  it  has  been  inserted  into  the  breech  of  a  "hot"  tank  gun.  The  temperature  of  the 
tank  gun  increases  as  consecutive  rounds  are  fired.  This  analysis  focuses  on  the  temperature 
profile  within  the  last  round,  after  it  is  inserted,  but  before  it  is  fired.  It  is  desirable  to  understand 
more  about  the  changes  to  the  propellant  temperature  under  these  circumstances,  as  these 
changes  may  affect  the  accuracy  of  targeting. 

MODEL  ASSUMPTIONS 

A  model  has  been  developed  to  represent  the  propellant  temperature  distribution  as  a 
function  of  time.  The  special  features  of  the  model  are  as  follows: 

•  The  chamber  and  the  round  are  considered  to  be  "long,"  so  that  the  conduction  is 
purely  radial  through  the  breech  wall  and  through  the  round  charge. 

•  The  heat  transfer  on  the  outer  surface  of  the  breech  is  due  to  free  convection. 

•  The  ambient  is  the  interior  environment  of  the  turret. 

•  The  charge  temperature  is  uniform  throughout,  upon  insertion  into  the  breech. 

•  The  thermal  properties  of  the  charge  are  the  best  available  mathematical  estimate  for 
particulate  charge  materials  of  the  M829A2  ammunition,  which  uses  JA2  propellant. 

•  The  initial  radial  temperature  distribution  in  the  chamber  was  based  on  the  Benet 
Technology  Team  data,  which  also  allowed  for  axial  variations.  The  "worst  case" 
axial  location  was  chosen  for  our  model. 

•  The  thermal  resistance  of  the  gap  between  the  round  and  the  barrel  was  not  modeled. 
The  results  represent  the  upper  bound  prediction  of  the  heat  going  into  the  propellant 
bed. 

MATERIAL  PROPERTIES 

The  model  contains  a  number  of  variables  that  represent  properties  specific  to  JA2 
propellant  and  the  properties  specific  to  gun  steel.  The  charge  and  steel  variables  of  interest  are 
density,  specific  heat,  and  thermal  conductivity.  The  values  for  gun  steel  were  readily  available 
at  Benet;  however,  finding  information  regarding  JA2  proved  to  be  much  more  challenging. 

After  much  research,  values  have  been  obtained  for  solid  JA2  material.  The  JA2  in  this 
ammunition  is  in  granular  and/or  stick  configurations.  The  values  for  solid  JA2  needed  to  be 
adjusted  mathematically  to  compensate  for  the  amount  of  air  (or  void  fraction)  contained  among 
the  small  pieces  of  the  propellant.  The  void  fraction  was  estimated  to  be  10%  in  this  analysis. 

The  solid  values  of  density  and  specific  heat  have  been  modified  using  this  value.  The  equations 
for  the  void  fraction  analysis  are  as  follows: 
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fa  =  vyv, 


(i) 


The  mass  of  the  charge  is 


tnc=PaVa+pcVc  (2) 

Using  the  above  two  equations,  the  average  density  of  the  charge  is  derived  as 

Pc=Pafa+Pc(l-fa)  (3) 

The  average  specific  heat  of  the  charge  is  derived  in  a  similar  manner  and  is 

ZP  =  CVoPafa  +  CpcPc(  l~fa)/(pja+(l- fa  )Pc  )  (4) 

The  mean  diffusivity  is  then 


ac=kc/(pcCpc) 


(5) 


where 


Pccp  C  paPaf a  +  C  pcPc(  1  fa)  (6) 

During  a  previous  project  at  Benet  Laboratories,  a  value  for  thermal  conductivity  for 

granular  JA2  was  obtained  experimentally  as  0.0774  Btu/(hr-ft-°F).  This  value  was  used  in  this 
analysis. 

MATHEMATICAL  MODEL 

As  stated  previously,  the  mathematical  model  of  the  propellant  is  one-  dimensional, 
where  the  radial  position  is  the  independent  variable.  The  model  uses  implicit  finite-difference 
equations  and  a  fixed  time  step.  The  propellant/gun  radial  distance  is  broken  into  a  number  of 
nodal  positions  as  seen  in  Figure  1.  The  nodes  are  characterized  by  equations,  which  can  be 
categorized  as  follows: 

•  The  center  node 

•  The  interface  node 

•  The  exterior  node 

•  The  interior  nodes 

The  center  node  is  the  radial  center  of  the  charge.  The  interface  node  is  the  radial 
location  of  the  interface  between  the  round  and  the  gun.  The  exterior  node  is  the  radial  location 
of  the  outside  diameter  of  the  breech.  The  interior  nodes  are  all  other  radial  nodes  located  either 
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within  the  propellant  or  the  gun.  The  nodes  are  equally  spaced  by  A r  with  one  exception;  the 
distance  from  the  center  node  to  the  next  node  is  Ar/2. 

Using  conservation  of  energy,  the  heat  going  into  each  node  is 


Qin  —  Q storage^  Q  out  (7) 

The  equations  for  the  various  nodes  are  derived  by  applying  the  conservation  equation  to 
each  elemental  shell.  The  derived  equation  for  the  interior  shell  nodes  is 


T’  =(l+2Fox  W:" -Fo,lU*(Arm  /2r.  W”  +<!-<*. /2r.  )jr£  )  (8) 


where 


Fox=(axAt)/ArJ  (9) 

and  x  =  a  variable  that  represents  either  "b"  for  the  gun  steel  nodal  region  or  "c"  for  the 
propellant  nodal  region. 

The  equation  for  the  exterior  atmospheric  boundary  node  is 

T:  +2Fo„BiT,  =lH-2FotBi  +  2Fot(l-(Ar„/2r„  ))]T 
-2FoO-(Arm/2r.))T£ 

where 

Fob=(ccbAt)/Arm2  (11) 

and 

Bi=(hArm)/kb  (12) 

The  equation  for  the  interface  shell  node  is 

(l/2)[(l/Foc  )+(k„/kc  )(l/Fob )lT'=t(l/2)l(l/Foc )+(kb/k,  ) 

(l/Fo„  )]+[<l-(Ar„ /2rm  ))+(k„  /k,  Kl+(Ar„ /2rm  BUT'"  (13) 

-(l-(Arm /2r.  )jr£  -(k„ /k,XI+(Arm / 2r. 

The  equation  for  the  center  node  in  the  charge  is 

Tl,  =(l+4Fo ,  wrl -4FOJT'  (14) 
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The  number  of  equations  is  equal  to  the  number  of  nodes.  The  equations  may  be 
reorganized  and  put  in  matrix  row  form.  The  equations  would  then  form  the  matrix  equation  as 

[AJ[TJp+i  =  [T]p  (15) 

Inverting  the  A  matrix  and  multiplying  through  yields 

[T]p+I  =[A]~1[T]p  (16) 

This  formula  is  used  iteratively  in  the  thermal  emulator  software  program  to  compute  the 
radial  temperature  distribution  in  the  charge/gun  model  as  a  function  of  time.  The  results  of  this 
program  were  validated  using  the  Heisler  charts. 

The  initial  temperature  distribution  in  the  breech  wall  is  an  important  feature,  which 
affects  the  accuracy  of  the  model.  This  distribution  depends,  in  part,  upon  the  firing  scenario 
prior  to  insertion  of  the  round  and  the  ambient  temperature  conditions.  Benet  has  a  software 
program,  FDHEAT,  that  profiles  the  gun  temperature  distribution  when  given  the  type  of  ammo, 
the  type  of  gun,  the  firing  rate,  and  the  ambient  temperature  conditions.  The  results  of  this 
software  program  are  used  as  initial  conditions  in  the  input  temperature  matrix  of  the  thermal 
emulator  program. 

The  average  temperature  of  the  propellant  with  respect  to  time  is  be  derived  from  the 
temperature  distribution  as  it  varies  with  time.  The  equation  for  the  energy  stored  in  the 
propellant  is 


ET=AE,  +  "f^AEm+JEl 

m=c+2 


(17) 


The  average  temperature  of  the  round  at  any  point  in  time  is  be  calculated  as 

Tav=Ti+ET/(pccpcVt)  (18) 

The  average  temperature  of  the  charge  as  a  function  of  time  may  be  useful  to  the  fire  control 
engineers. 
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RESULTS 


Figures  2  through  6  present  a  sample  of  runs  including  propellant  average  temperature 
versus  time  for  various  initial  round  temperature  and  gun  temperature  scenarios.  Each  graph  has 
a  plot  for  the  average  temperature  in  the  last-loaded  round  after  a  certain  number  of  rounds  is 
fired.  All  scenarios  have  a  constant  rate  of  fire  of  six  rounds  per  minute. 

The  graphs  indicate  that  a  heated  chamber  definitely  causes  a  change  in  the  average 
temperature  of  the  propellant.  The  heat  transfer  is  significantly  affected  by  a  number  of 
variables,  including  the  number  of  rounds  previously  fired  and  temperature  of  the  round  obtained 
from  stowage.  Figure  6  indicates  that  the  convective  heat  transfer  coefficient  plays  a  significant 
role  in  that  circulating  air,  as  opposed  to  stagnant  air,  within  the  crew  compartment  increases  the 
cooling  process  of  the  gun  and  therefore  affects  the  propellant  temperature. 
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Constant  Axial  Location 
Rate  of  Fire  =  6  rounds/min. 

Ambient  Temperature  =  70°F 
Propellant  Initial  Temperature  =  35*F 
Number  of  Rounds  Fired:  6, 12, 18, 24, 30, 36 
1.0 


Rds 


Figure  2.  Propellant  average  temperature  versus  time  with  mild 
crew  compartment  temperature  and  cold  round  temperature. 
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Average  Temperature  of  Propellant  (°F) 


Constant  Axial  Location 
Rate  of  Fire  =  6  rounds/min. 

Ambient  Temperature  =  70°F 
Propellant  Initial  Temperature  =  70°F 
Number  of  Rounds  Fired:  6, 12, 18. 24, 30, 36 
h-  1.0 
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Figure  3.  Propellant  average  temperature  versus  time  with  mild 
crew  compartment  temperature  and  mild  round  temperature. 
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Average  Temperature  of  Propellant  (°F) 


Constant  Axial  Location 
Rate  of  Fire  =  6  rounds/min. 

Ambient  Temperature  =  70°F 
Propellant  Initial  Temperature  =  90°F 
Number  of  Rounds  Fired:  6, 12, 18, 24,  30,  36 
h  =1.0 


Figure  4.  Propellant  average  temperature  versus  time  with  mild 
crew  compartment  temperature  and  hot  round  temperature. 
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Average  Temperature  of  Propellant  (°F) 


Constant  Axial  Location 
Rate  of  Fire  =  6  rounds/min. 

Ambient  Temperature  =  90T 
Propellant  Initial  Temperature  =  90*F 
Number  of  Rounds  Fired:  6, 12, 18, 24,  30,  36 
ft  =  1.0 


Figure  5.  Propellant  average  temperature  versus  time  with  hot 
crew  compartment  temperature  and  hot  round  temperature. 
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Average  Temperature  of  Propellant  (°F) 


Constant  Axial  Location 
Rate  of  Fire  =  6  rounds/min. 

Ambient  Temperature  =  70PF 
Propeliant  Initial  Temperature  =  70*F 
Number  of  Rounds  Fired:  6, 12,  18,  24,  30,  36 
h  =1.0  and  h  =  5.0 

36  Rds 
30  Rds 
24  Rds 

18  Rds 
12  Rds 


6  Rds 


*  *  *  *  '  '  ’  1  ^Notc:  Solid  lines  h*=  1.0 

0  1800  3600  5400  7200  9000  10800  12600  14400  Broken  lines  ft- 5.0 

Time  Since  Insertion  of  Round  (sec) 

Figure  6.  Propellant  average  temperature  versus  time  with  mild  crew 
compartment  temperature  and  mild  round  temperature  at  h  equals  1.0  and  5.0. 
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NOMENCLATURE 


The  variables  for  the  model  are  assigned  as  follows: 
m  Node  number 

p  Beginning  of  time  step 

p+1  End  of  time  step 

A rm  Radial  distance  between  consecutive  nodes 

r  Radial  distance  starting  from  center  of  charge 
T  Temperature 

[T\  Temperature  matrix 

mc  Mass  of  charge 

kc  Thermal  conductivity  of  charge 

h  Free  convective  heat  transfer  coefficient 
h  Thermal  conductivity  of  breech  (gun  steel) 

Bi  Biot  number 

a  b  Thermal  diffusivity  of  breech  (gun  steel) 

Fox  Fourier  number  of  material  x 
a  c  Thermal  diffusivity  of  solid  charge 
ac  Average  thermal  diffusivity  of  charge 
cpa  Specific  heat  of  air  in  the  charge 
cpc  Specific  heat  of  solid  charge 
cpc  Average  specific  heat  of  charge  -  air  mixture 
p  a  Density  of  air  in  the  charge 
p  c  Density  of  charge 
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pc  Average  density  of  charge  -  air  mixture 

A  t  Time  step 

q  Heat  flux  rate 

V,  Volume  of  solid  charge 

Va  Volume  of  air  in  charge  -  air  mixture 

Vc  Volume  of  solid  JA2  material 

fa  Void  fraction  of  charge 

Ec  Energy  contained  in  the  center  shell  of  charge 
Ei  Energy  contained  in  the  outer  shell  of  charge 
Et  Energy  stored  in  the  charge 
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